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CHAPTER  1 


INTRODUCTION 

Any  material  which  has  three  mutually  perpendicular 
planes  of  elastic  symmetry  is  said  to  be  orthotropic.  If 
from  such  a material  a thin  flat  rectangular  plate  of  uni- 
form thickness  is  cut  such  that  the  plane  of  the  plate 
is  parallel  to  a plane  of  elastic  symmetry,  then  this  plate 
will  have  two  perpendicular  axes  of  symmetry  in  the  plane 
of  the  plate.  Such  a plate  is  said  to  be  •orthotropic.  In 
this  dissertation  the  state  of  stress  in  a large  rectangu- 
lar plate  of  this  type  containing  a circular  hole  filled  by 
an  Isotropic  material  is  considered.  The  problem  is  an 
important  one.  As  an  Illustration,  the  problem  of  a wood 
plate  containing  a circular  knot  can  be  cited  since  both 
plain-sawn  or  quarter-sawn  boards^  are  examples  of  orthor 
tropic  plates.  In  the  discussion  the  large  olate  will  be 
assumed  to  be  subjected  to  a uniform  tension  S at  large 
distance  from  the  isotropic  disk  and  directed  parallel  to 

^Alvin  Jev/el  Owens,  "Effect  of  a Rigid  Elliptic  Disk 
on  the  Stress  Distribution  in  an  Orthotropic  Plate." 
Doctor's  thesis.  University  of  Florida,  1950,  p.  2. 
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an  axis  of  elastic  symmetry.  The  stress  functions  for  the 
orthotropic  plate  and  for  the  isotropic  disk  will  be  ob- 
tained and  an  examination  will  be  made  of  the  stresses  near 
the  disk.  It  is  assiomed  that  the  strains  always  remain  with- 
in the  limits  of  perfect  elasticity,  l.e.  the  material  always 
returns  to  its  initial  shape  after  the  external  loading  is 
removed. 

For  the  orthotropic  plate  mentioned  above  with  no  holes 
and  with  a uniform  tension  applied  along  tv/o  opposite  edges 
there  are  no  stress  concentrations.  This  means  that  the 
ratio  of  the  normal  stress  at  any  point  in  the  direction  of 
the  applied  tension  to  the  normal  stress  at  infinity  is  one; 
the  ratio  of  the  normal  stress  at  any  point  in  the  direction 
perpendicular  to  the  applied  tension  at  infinity  is  zero; 
and  the  ratio  of  the  shear  stresses  at  any  point  in  either 
of  the  above  directions  to  the  normal  stress  at  infinity 
equals  zero.  As  has  been  stated  the  values  of  the  ratios 
given  above  are  constant  for  all  points  of  the  plate  so  that 
there  are  no  stress  cone 'ntrations  anywhere. 

When  a disk  is  Inserted  in  a hole  in  an  infinite  ortho- 
tropic  plate  and  a uniform  tension  is  applied  along  opposite 
edges  of  the  plate  as  in  the  present  case,  then  the  ratios  of 
the  stresses  at  a point  to  the  normal  stress  at  infinity  vary 
from  point  to  point.  There  are  definite  stress  concentrations 
in  the  neighborhood  of  the  disk,  that  is,  the  ratios  of  the 
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stresses  at  points  in  the  vicinity  of  the  disk  to  the 

normal  stress  at  infinity  are  larger  than  the  corresponding 
ratios  for  some  other  points  in  the  plane.  An  examination 
of  the  stresses  in  the  vicinity  of  the  disk  will  be  made 
with  tne  idea  in  mind  of  determining  the  points  of  greatest 
stress  concentrations. 

Dr.  A.  J.  0wens2  treated  the  problem  of  an  infinite 
ort'ootrooic  plate  with  a small  elliptic  h61e  at  the  center 
filled  With  a rigid  elliptic  disk.  In  this  problem  a circu- 
lar hole  In  an  infinite  orthotropic  plate  is  filled  with  a 
circular,  non-rigid  disk  of  an  isotropic  material. 

Professor  C.  B.  Smith-^  treated  the  problem  of  an  infin- 
ite rectangular  orthotropic  plate  with  a small  elliptic  or 
circular  hole  at  its  center.  Professor  Smith  found  greater 
stress  concentrations  near  the  free  hole  than  Dr.  Owens  did 
near  the  rigid  elliptic  disk.  For  this  problem  the  stress 
concentrations  are  found  to  be  intermediate  between  those 
found  by  Professor  Smith  and  those  found  by  Dr.  Owens. 


^Owens,  OP.  clt. 


^C.  n.  Smith,  Effect  of  Elliptic 
the  Stress  Distrioutlon  in  Plates  of  V 
ered  asOrthotropic  Materials. 

Report  No.  iblO,  194i+.  " 


or  Circular  Holes  on 
Vood  or  Plywood  Consid- 
Products  Laboratory 


CHAPTER  II 


SOLUTION  FOR  PLANE  STRESS 


The  center  of  the  circular  disk  of  radius  a Is  chosen 
as  the  origin  of  the  coordinate  system.  It  is  assumed  that 
the  hole  Is  at  the  center  of  the  plate (whenever  the  word 
plate  Is  used  the  region  exterior  to  the  disk  is  intended). 
The  axes  are  taken  to  lie  parallel  to  the  two  axes  of  elastic 
symmetry  of  the  plate,  hence  the  uniform  tension  at  infinity 
acts  parallel  to  a coordinate  axis.  (Pig.  1) 

The  mathematical  notation  of  Professor  S.  Timoshenko^ 
v;ill  be  used  throughout  this  treatise.  Using  this  notation 
the  boundary  conditions  on  the  periphery  of  the  plate  are: 


<1^  ( X = ± oo)  _ 5 


(Y  - ±»)  - O 


(1) 


ii 

Timoshenko,  S.,  Theory  of  Elasticity  (New  York:  McGraw- 
Hill  Book  Company,  1951 ) • 


hr 
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The  boundary  conditions  on  the  boundary  of  the  disk  are, 
using  a prime  to  denote  the  stresses  and  displacements  in  the 
disk: 


The  plate  is  considered  to  be  thin  and  the  thickness 
in  Z direction  is  taken  equal  to  unity.  Since  the  X-  and 
Y-  axes  are  the  axes  of  elastic  symmetry  of  the  plate  and 
since  the  plate  is  in  a state  of  plane  stress  the  follow- 
ing relations  between  the  components  of  stress  and  strain 


(2) 


► n = a 


u - u 


V'  V 


hold: 


dU 
' dX 


X 


(3) 


In  these  equations  and  are  Young’s  moduli  in 
the  X and  the  Y directions  respectively.  The  symbol  is 
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Poisson's  ratio  and  xj  ratio  of  the  contraction 


parallel  to  the  Y-axls  to  the  extension  parallel  to  the  X- 
axls  associated  with  a tension  parallel  to  the  X-axis.  The 
symbol  G is  the  modulus  of  rigidity;  G^y  is  the  modulus  of 
rigidity  associated  with  the  directions  of  X and  Y, 

For  plane  stress  the  equations  of  equilibrium  are: 


These  equations  of  equilibriurii  are  obtained  from  the 
conditions  of  static  equilibrium  without  reference  to  the 
laws  relating  stresses  and  strains.  They  obviously  will  be 
satisfied  by  any  stress  function  F such  that 


d Txy  , ao^  - 

3X  ay 


2. 


(5) 
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The  compatibility  equation  for  plane  stress 

3 Sx  , d ^ 

dY^  " 3X3Y 

is  obtained  by  utilizing  the  notion  that  u(x,y)  and  v(x,y), 
the  displacements  in  the  X and  Y directions  respectlvely^are 
single-valued  and  continuous. 

Substitute  equations  (5)  into  equations  (3)  and  obtain 


11 

X 

(L) 

1 dV 
E,  3Y" 

v*yx  a F 
Ey  SX^ 

CD 

1 1 

1 d^F 

Ely 

U«Y  d^^F 

(7) 

II 

>- 

-1  d"F 
axdY  . 

Using  Vyx 

Vxy  equations  (7)  become 

" Ex 

1 

VxY 

E,  9Y' 

3X^ 

Y 

1 aT 

Ey  3Y"- 

(8) 

dXdY  , 
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Substituting  these  values  into  the  compatibility  equa- 
tion, there  results 


_L 


^ -V 


4- 


I 


■r 


dV_v;  dT 


- -I  dV 
G;,YdXW 


ax'* 


-I- 


_L 

G 


XY 


2Uxr 

E. 


aT  , 

BX^dY^  * 


E,  a-f 


(9) 


Equation  (9)  can  be  written  as 


dV 

dX*t 


+2K|£;.-|-  ■?¥.  = O 


ax'ari' 

c; 

where"^  I 

K-  ttv 

T\  r gY 


arv'*' 

\ 2.  VxY 


G 


xY 


E- 


and 


WE> 


(10) 


Substitution  of  P “ P(X■^S'^_  ) Into  the  new  form  of  the 
compatibility  equation  (10)  leads  to  the  following  equation 
in  5 

8'^-)-2KS’'+i  ~ O. 


^K.  W,  March,  Flat  Plates  of  Plywood  under  Uniform  or 
Concentrated  Loads,  Eorest  Products  Laboratory  Report  Uo.  1312, 
Marcn,  19 i+2,  p.  q.0. 


10 


This  is  a quadratic  in  that 

Since  for  wood  K defined  is  probably  always  greater 
than  one^,  let  ~ cosh  ij)  then 

— cosh  (j)  ± J cosh^ (J)  — • I 


S'r 


cosh  (()  ± sinh  (j) 


Therefore  the  roots  of  equation  (11)  are 

S,~  te  ^ 

ce  ^ 

S3  = -ie"^ 

6 = . 

Uv|,  — , • 


The  compatibility  equation  (10)  can  also  be  rewritten 


as 


dX- 


3^  , 


3^1  3^2. 


F=o 


wnere 


2 

OC  IT  C, 


; 


e 


-<t> 


s _nce 


e^4-  z: 


2.  cosh  i'  — z K . 


6 


Smith,  cit.  , p. 
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The  differential  equation  (12)  is  satisfied  by  any 
stress  function  of  the  form 


The  stress  function  in  (13)  may  easily  be  demonstrated  to  be 
a solution  of  (12)  by  direct  substitution. 

Thus  equations  (12)  and  (13)  are  valid  for  the  stress 
function  associated  with  any  problem  of  plane  stress  in  an 
orthotropic  plate  with  the  X-  and  Y-axes  being  axes  of  elas- 
tic symmetry. 

For  the  problem  \mder  consideration  a, stress  function 
of  the  form  (13)  which  satisfies  the  boundary  conditions  at 
infinity  and  on  the  boundary  of  the  disk  is  sought. 

7 

Professor  Smith  found  in  his  problem  of  an  infinite 
orthotropic  plate  with  an  elliptic  or  circular  hole  at  the 
center  that  a satisfactory  stress  fimctlon  is 


(13) 


where 


^ ^ (z.+ w.)] 


2. 


Smith,  0£.  clt. 
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where 

z,  - x + 


Z-^=  x + i-^eY 

W2.=/zF^ 
a-^Vb^  . 


It  is  recalled  that  H ;;eY 


the  hole  is  given  by 

X - a cos  © 
Y r b sin  © 


Also  the  boundary  of 


and  a uniform  tension  is  apolied  at  the  ends  of  the  plate. 

Owens^  used  this  same  stress  function  for  his  problem 
which  has  been  mentioned  previously.  Since  Owens'  stress 
function  satisfied  the  boundary  conditions  at  infinity  it 
seems  likely  that  his  stress  function  will  be  suitable  for 
this  problem  if  it  satisfies  the  boundary  conditions  on  the 
periphery  of  the  disk.  It  is  found  that  his  stress  function 
does  satisfy  the  imposed  boundary  conditions. 

The  stress  function  then  is 


(lii) 


Qowens,  0£.  cit. , p.  l5. 
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where  A = A^^-1-  IA2 


B = + IB2 


Z|=X+i.=ieV  + 


but  where  now 


ere  now 

_ q'h'i  oHl-/5-e^)  ■ 


The  boundary  of  the  disk  is  given  by 

X n a cose 
Y* r Q sine 

and  S is  the  uniform  tension  applied  to  the  ends  of  the 
plate  in  the  X direction. 

To  insure  that  the  stresses  be  single-valued  the  square 
roots  in  (If?)  so  chosen  that  the  inequalities 


always  hold.  After  the  stresses  are  obtained  it  will  be 
more  evident  why  such  restrictions  are  necessary.  There- 
fore after  the  expressions  for  the  stresses  are  obtained 
these  inequalities  •will  be  discussed  further. 

For  convenience  of  reference  a list  of  useful  relations 


will  be  given  here. 


Ik 


- af. 

az. 

_ af, 

ax 

~ B'2-. 

dX 

B2, 

dfzCA) 

^ afi. 

d Z-^ 

_ dh 

dX 

" aZi 

ax 

az^ 

_ df, 

az>, 

BY 

~ 5Z, 

3Y 

_ Bfi 

aZx 

- afi 

■ay 

- 

SY 

3^2 

aw/  __ 

H, 

3 2,  " 

W, 

"aWi. 

3ZA 

'2. 

rL. 

Wj,  , 

F’rom  the  stress  function, using  the  above  relations  where 
necessary,  one  has 


dX  '' 


A 

~B 

ZTsi 


fe-w:)0- 


(zr^2)(i- 


W2- 


Rev;rite  as  follows 


a£=R 

ax 


Since 


then 


A_ 


VV,(E,^-vV)  W.(z,tw,) 


B 


w,  = 7z,‘-r,^ 
zr-  w,"  = . 

BF 


Using  this  relation  -Z-h  becomes 

3 X 


9F  _ 


ax 

or 

aF 

ax 


R 

A 

'-cz-.-w;)  ■ 

, B 

-(Z^-W^)  + fz^+Wz) 

r\ 

z 

W,  ( Z,  t W.) 

Wz(H.x  + Wx) 

-R 


_ sV 


A 


Z-.-t-W, 


~4~ 


B 


3 X 


-R 


A 


Zj-fvV^ 

-C+lr)  ■ 


^ B 


- ( I + \A/.  ^ 


(Zi+w».V 


But  this  gives 


a;-  R 


-A 


-B 


W.  ^Z,+-W,) 


Wi(z.  + w^') 


(i6) 


It  now  becomes  apparent  what  the  form  of  the  express- 
ion for  the  other  stresses  will  be.  Performance  of  the 
necessary  differentiation  yields. 


o;  = 


a^F 

3Y 


^delA_  ^ 


WYz,+w,)  W^(Z^t-W^) 


+ s 


(17) 


16 


and 


3X9Y' [ W,(£^+W!)  Wi(?^+VVi,) 


(18) 


The  stress  formulas  (l6),  (17)»  and  (l8)  contain 
multiple-valued  expressions 


The  branch  points  of  the  function  W are  readily  seen 
to  be  ^ 1C.  Since  there  are  two  and  only  two  branch  points 
the  branch  cut  may  be  taken  as  a straight  line  joining  the 
points  1C  and  -1C.  How  does  this  branch  cut  map  onto  the 
W-plane? 


w.fz,  +w.) 


Wz(Z2.+  w'i)  ’ 


Assume  that  ^ O 
and  1st  -~C^, 


Define 


w = z.+w,  r z,  ± / zt + . 


Given  H, 

X.,  = X-O 


and 


along  the  branch  cut 

W r u + 
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where 

whereby 


U = ± y C^-  Y,^ ' ; V = Y,  . 

^ m a circle 

with  center  at  the  origin  in  the  VV  plane  and  of  raduis  C-| 
If  we  restrict  W such  that 


u"  + 


|W|  = |Z,+ W,|  >y\c\-\^> 


then  the  variable  % is  prohibited  from  following  a path 

which  crosses  the  branch  cut  and  thus  W remains  single- 
valued. The  entire  plane  with  the  exception  of  the 

branch  cut  is  mapped  into  the  region  in  the  \H  plane  exterior 


to  the  circle  at  the  origin  of  radius 
Assume  > O then 


c 


Now  the  branch  points  of  IV  are  db.^1  ^ and  a logical  branch 
cut  is  a straight  line  joining  the  points  iT,  and  - "iT  , in 
the  2.,“  plane.  The  branch  cut  maps  into  a circle  in  the 
VV-  plane  with  center  at  the  origin  and  of  radius  . The 

restriction  used  before 


|wl=iz.  + w, 


Insures  that  the  branch  cut  will  not  be  crossed  and  that  W 
will  remain  single-valued. 

Obviously  an  entirely  analogous  argument  relative  to 
and  yields  the  relation 
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+ W, 


in  order  that  the  stresses  remain  single-valued. 

It  is  necessary  to  show  that  the  given  stress  function 
yields  stresses  v;hich  satisfy  the  boundary  conditions  at 
infinity.  When  X - ±.  then  "Z,— 

Z.+w,  - • 

Examining  equations  (17)»  and  (l8),  it  is  found  that 


o; 


X-> 


,]=R[ 


O -t-O 


+ s = s 


provided  that  the  constants  A and  B are  not  infinite. 
Likewise 

- o. 

For  the  two  edges  of  the  plate  Y = 

and  ^,+w.  = Zi+Wz  = ^ • 

Equations  (l6)  and  (l8)  yield 


T 

' X Y 


X -=►  - OO 


= R 


0 + 0 


Y->  CD 
Y->-  oo 


Y-^  oo 
Y ^-OD 

m 


R 

R 


o + o] 

o + ol 


r O 
= O . 


Hence  it  is  seen  that  all  the  boundary  conditions  at 
the  edges  of  the  plate  are  satisfied. 
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There  remain  then  four  boundary  conditions  to  be  satis* 
fied  at  the  periphery  of  the  disk.  These  four  conditions 
mathematically  stated  are: 


Tn©  n Trt( 


n-Q 


(2) 


u - u 

v'  - V 

where  the  primes  denote  the  disk. 

Since  can  be  obtained  from  the  rectangular  stresses 
CT^  , and  Txv  the  relation  between  them  v/111  be 
obtained. 


Figure  2 

C.oS^©  -t-  cr^  5'n*0  -f-  2 s>'n©  cos©  . 


(19) 
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Substitution  of  equations  (16),  (1?)  and  (l8)  into 
equation  (19)  yields 

+ — ^ Jj3^€\os'g  -sin*e  + 5i.n0  cos0^ 

nJx(2z+vVi)'^ 


-H  S cos^e  . 


On  the  boundary  of  the  disk  r = a.  »Vhat  is  the  value 
of  and  V/^  when  r - a? 

Since  X - a cos  9 and  y = a sin  9 then 

2, +W,  — C (co5o  ^ L c(e  5ifi e) 

^ Ja^(co5^o  * 2i«Ke5/>?o  cose-ct^^^-^io^^aYz-at^e^) 

or  Z.  + W.-  orco50^-^*ote  5/n0) 

/_  5in'd  + cose  - o(^^sin^e  ^ . 


Using  a trigonometric  identity  this  becomes 


+ 


VV,  — C?  ( C05© -h  < S/O  0) 


^ 0 +-  ^'*^0  C03  0 + o<^€\os*'€ 

2 f W - 0 ^coso  + t’oce  5in0  +^(<>(6  cose  + l 5I  ne)^ 
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Thus 


or  rewriting 


-l-W|  — O ( I ^ ci6)( CoS  e -h  <■  5 i R , 


Prora  the  above  computation  it  is  seen  that 


W,  =:  a (ocC  coSQ  ^ L sine)  ^ 

Similarly 

W^ii  a(Se  cose  f iSin  e) 
and  a(\  + se)(cos  e + isin^) . 

The  above  values  for  VV2»  +■  , and  &2  ''^Z 

permit  the  rewriting  of  C3^  . 


crs(n  = <a)  r R 


(oc€  C05©  -e  t 


+ 5 CCX5*e  . 


This  reduces  to 

A(o^6co5e  +j  smeXcose-tscyic) 
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Inasmuch  as  It  will  be  desirable  to  have  expressed 

in  terras  of  the  sines  and  cosines  of  various  multiples  of 
the  angle  0,  the  powers  and  products  of  the  sine  and  cosine 
functions  will  be  converted. 

(o<'e  Cos^  f i (l  -o(fe)  5in  © cos  © + ^ i f\j^ 

e co^^e  + ici-'Se) sine  cose  +sin^e 

a V/  f 


) 


a\(n^a) 


-I- 


5 


© . 


By  actually  taking  the  real  part  of  the  above  expression 
there  results 


^,(o6fcco3 qh- stn^e) - {\r.(t~oLe)  sine  cos  © 

C / v-ote) 


B,  (^6.  cos^o  tSLn^o)  + 3u(Be-i)si no  C.OSQ  i 

+ — ^ ^ CoS  e . 

a^C  I +3€) 


Now 


Stn©  Cos©  - 


S ■-  n 2.Q 

Z 


CoS  © - 


+ cos  2.0 
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, __  1-COSZ0 

and  S I r\  © ~ — 

Z 

These  relations  perinlt  the  expression  of  CT^  in  the  fol- 
lowing form: 

Q-  - ^Xo<e(^^^co5  ze  i)-  cosze]  P\A‘>(e-\)  stn  ze 

Za^(  i-i-<=(e)  2q^(ii-^€:) 

cos  •z.eH  /-  coszj  B^(s^-()sitiZ9 ^ ^ S cosze 

Za'*‘(  i-hSe)  %a^(nse)  ^ ^ 2. 


This  can  be  rewritten  as 


^ Z^V/y-o<6; 


_! ^^/•/•elre)A,  cosze  +<<<e-i)f\^s.n2« 

V/y- 0(6)1-  J 


. (^6-/)6,co52©-4-('^6-l)e^5/n2©J 

Z^7Y/f^€)L 


(20) 


, S . S cosze 
^ ~Z  ^ 2 


Next  an  expression  for  the  shearing  stress  Tn»  in 
the  above  form  v/ill  be  obtained  from  the  rectangular  stress 
components.  Utilizing  Figure  3 it  is  readily  seen  that 


\e  ~ ~cr^  Sin  © cos  e -i-  Oy  sine  cose  + - sin*©)  . 


Piglire  3 

Use  expressions  (l6),  (17),  and  (18)  to  obtain 


f si'ngcose-nVeL^iVe 


a^(i  -t  o(eX°<^  cos  o + t 5 4 ne) 


_ S{cos  0 - 1 Sen  ^ Q t Stne  cose  ■j-i^^^^'Oo-coso)) 

CJ  ( I C0SO+  c Sirie) 


~~  S si  r\&  <cose  . 


Tnls  will  slrnpllfy  somewhat  by  factoring  to  the  form 


%0=  - 


R 


/\(cos 0-4.510  0^(^6  Cose  »*'S<ne)(^65<n^- 


e COSO  i- 1.  si.r\0 


a^(l-tc<€) 
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■w 


l 

R 

B Ccoi  e-  L 5 ( a G)C/«e  cos«  4 isln^Of  -/r«< 

a^(/-tse)  ' ' 

_ Ssinxe 

~z 

eecosQ  -f  t’stn  e 

- 

--kR 

(JK\\  1 h; 

Ly°<€Smeco5e-il*€Stn^e  + cose]-  smo 

'ne  (x^ 

1 + o<€ 

* fetiB0(^65inOcO5e-t(^65iV©4  Co5*-0)  — 5in&  Co<.nS] 

S Sin20 

liae 

Z 

Now  actually  take  the  real  part  of  the  expression,  and  one  has 


-(°<^-<)stnecos0A.  -^o<£Sine->-coS^Q)Ra- 


^ Sinoco^e  B,  — (g€5/ng  ■»  cos^e)  Ba 

0'(m  S€) 


S sinze 


Apply  the  trlgonoraetric  identities  previously  used  in 
obtaining  \ it  follows  that 


"7  - 

‘n©  - 


n 26  — [ I 4- oce +- (>-cce)cos  i©] 
Zo^  C I +ot60 


(21) 

_ (56 -/)b,s in  2e-t-[;u^£4(r-4g6)cosze]B2  Ssinz.© 
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Now  the  polar  stresses  and  ^riO  for  the  plate 

are  obtained  in  the  desired  form.  In  order  to  satisfy  the 
boundary  conditions  at  the  periphery  of  the  disk  it  is 
necessary  to  have  the  corresponding  stresses  for  the  disk 
and  for  these  stresses  to  be  in  the  form  indicated  above. 
V/hat  is  the  stress  function  for  the  disk? 

S.  Timoshenko^  has  given  the  general  solution  for  the 
two  dimensional  problem  of  an  isotropic  ring  or  portion  of 
a ring,  which  includes  as  a special  case  the  problem  of  the 
disk.  This  solution  is  given  in  the  form  of  a stress  func- 
tion and  is  set  forth  below  with  the  prime  indicating  that 
the  expression  refers  to  the  disk. 


fr'-  ajnn  + t Coff-I nn  -t-  don^o  + 9 

-f^nesine  -t- H b,  n In  ® 

_ Onecose  +-(<j,nhc;R-'tJ,'n  In  n) 

2 


OP  n 


-n  + 2. 


) COS  no 


n + z 

+• 


-nt2.\ 

n ) 5irt  ne . 


^Timoshenko,  op.  c it . , p.  lllj.. 
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To  prevent  laultiple-valued  stresses  take  - d - b 

O O O 

= dj  =0.  I'ere  the  primes  in  the  right  member  do  not  refe. 
to  the  disk. 

In  polar  coordinates  the  relations  between  the  polar 
stresses  (J and  tlie  stress  function  are 


n-' L ^ -L  £L 

^ - n.  an  A*-  3©' 


(22) 


1 I 


T- ' - ± ^ — i-  ^ F 

'ae  - n*'  *©  n anae  • 


2l£' 

’ ao  ’ ae' 


Thus  err  . ^ . ThlTL  , and 


arvae 


must  be  obtained. 


First  differentiate  P'  with  respect  to  A . 


= Z A t ^ Q,  A coSe  - 

^n.  2.  2 


+ (3d.n^-c!n'^)  ^«n0 


^ / n-i  , ,1  n + » ' -n  + l i'  -n  + i\ 

+ £(na„a  i-CAtZlb„A  -nCl^A  -Cn-Z^bnh  jeosn 


0° 

KT-  / ^ -'1''  , 1 r . i'  • 

+ +Cn  + Z)dnn  -nC^A  -{n-2)d„n  jsiAh© 


n=2. 


Since 


an 


is  multiplied  by --L—  in  the  expression  for 


and  since  r ranges  to  zero  then  it  is  obvious  that 
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a,,  a|,  c|  , a^,  c,!^,  and  drt  must  equal  zero  to  prevent 

infinite  stresses  at  the  origin. 

Thus  differentiating  with  respect  to  d and  noting  that 
the  above  constants  are  zero,  one  obtains 

I 

- C/o  - b,  n^6tr\©  - ^(co5  0- ©5i’n  e) 


- cl.^^CDSe  - sin  ne 

n=L  ' 


CO 


+ n C05  n © . 


}l~Z 


It  is  further  noted  that  Is  multiplied  by  in 


the  expression  for  ^0  . To  prevent  an  infinite  stress 
as  r — > o it  follows  that  a^  and  c,  must  be  zero. 

Next  — — - is  obtained  as  follows 


iE'  - 

ae^  - 


C05  ne 


nsz. 


00 


m-Z 


+ d^n"'")n"sin  ne, 

n-Z 
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One  other  partial  derivative  is  needed  and  is  obtained 
in  the  following  form 


^ / 

~ — 3 e -h  3 di  n^co5  0 


CO 


-E  (>ia„n"‘h  (n+2.)bnr\"*h  n sin  ne 
n=z  ’ 


-h  21  n cos  ne 


Substituting  these  values  into  equations  (22)  yields 


(Tf^-  2 bo  -»■  3 b, rt  Cos  0 i-  3 d,  n Sin  0 
-b^(  3- (n  + 2.)b„n”)  cos  ne 

n=l  ’ 

OP  j 

+E(ncnn"-  -H(n  + 2)  j 5in  ne  - b,n.cos  e 

n = Z 

— diftStnO  — XI  (^n  ^ + ^'cos  n© 

oo 

~Zr  (Cr^t\^  clryfi'^)  n^Si'n  n©  . 
n=z 
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A s Irapler  Torn  would  be 


CO 

n-'=  2h,,tZjp,nco5e  + Zcj,n^^ne  +2I(nq„n’^"+(n+z)kn'^co5n© 


CO 


(r\C^n^  i’  (n  + Z)drt^  )sin  ne 
n=z 


(23) 


o»  CP 

bnn")n\o5ne  ne  . 


n=7L 


n=z 


The  shearing  stress  for  the  disk  Is  found  to  be 

CD 

E 

n = z 


%'=  -h,n5in  e + cl,nco5e  — ZT  (apn"  +b„n"Jn  sin  ne 

/I  C7  ■ — ** 


OO 

)nco5  ne  +-  3b,risin©  -3cJ,;(  cob  © 

n = ^ 

^y(nQ^n'^ +(n+2.)bp^")n5<n  ne-ZK^^n^  )nco5ne 

JTz 


or 


Tn»  = Z-b.nsine - Zd,n  cose 


4-2^[(n-0Of,n"  Cn+0bn^^]n  s«’n  ne 

h=2. 

-f;.L(n-.)C„n"'i  (nn)d„n"]n  c.s  ne, 


(2lt-) 


nax 
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Observation  of  the  form  of  the  stresses,  (T^  , Xe. 
in  the  plate  reveals  that  these  expressions  contain  only 
constants,  and  sines  and  cosines  of  20.  Thus  from  equations 
(2)  the  stresses,  in  the  disk  must  contain  only 

constants  and  sines  and  cosines  of  29 o 

Comparison  of  equations  (20)  and  (23)  reveals  that  b, - 
d,  = 0 since  there  are  no  sin  0 or  cos  0 terms  in  CTl  , 


ft,  I B,  + ^ 

ZoVhoie)  za\n^Si)  Z. 


Therefore 


ua^  V 


(25) 


since  the  constants  of 
Comparison  of  the 


must  equal  constants  of  CTJ^ 
coefficients  of  cos  20  in  (20)  and 


(23)  yields 


_ (mie- 0 A.  , Ose-i)  B,  ^ _5_ 


rn 

X 


hus 


- _ (o^6-QA.  _ ^ 


(26) 
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Coraparlson  of  the  coefficients  of  sin  29  gives 

ZC-i_  + H - UC^-  ^ (d£-l)Bj^ 

ZQ^(sei-i) 


C — - Cs^-t)  B,  ^ 

^ — i/ij'^oce  + O 4^cj^C^e-n) 


The  same  procedure  is  followed  as  before.  A compari- 
son is  made  between  equations  (21)  and  (2i).).  Comparing  con- 
stant terms  indicates  that 


2a"  za'- 


or 


Likewise  comparison  of  the  coefficients  of  sin  29  gives 


-ZCz 


ZaYi+-<€)  ‘ 2a^(^?€+i)  z • 


But  the  expression  for  a^  is  known  so  that 

L ( / A>  __  ~ 0 _ 5_  . Ck6-0^«  ■ (se~i)B, 

‘4a^^'2aV«<€4  0 2^}N5e+i)  2 2aYo<«+i)  2Q‘(i5€ti;  2/ 


b,=  0. 
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Siinilarly,  equating  the  coefficients  of  cos  26.  yields 


-zC-^-U  a'di 


(i-o(€) Az.  _ 0-Se^^^ 


Substitute  here  the  value  of  Cj^^  obtained  previously,  then 


0-die)  __  (l-3e)  Pz. 


and 


A — (j — , (l 

iaHn-oLt)  4Q“(m-i5€)  ■ 


Thus  two  of  the  boundary  conditions  at  the  periphery 
of  the  disk  are  satisfied,  namely 

(Tf^ln  = a)  - cr;^  in- q) 

T„e(n^a)  = Xe(n^a). 

The  preceding  boundary  conditions  have  been  such  that 

^ I 

all  the  constants  of  p are  zero  except  b^  , a , c^  , 
and  dj  so  that  at  present  the  stress  function  for  the  disk 
is 

pr' -z  b„  Qin’’cos  2 0 + (c2^n^+ Stn  zo. 
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It  is  much  more  convenient  to  find  the  displac.ements 
In  the  plate  in  terras  of  the  rectangular  stresses  (J^  and 
0~y  so  the  rectangular  stresses  for  the  disk  will  be 
obtained.  This  can  be  done  very  readily  if  r and  9 are 
expressed  directly  in  terns  of  x and  y. 

Thus  r becomes 

r' = b*(xVY')  + Y"")  cos 

H-[c^(xVy^)  -h  d»(xVY^fjsin  2® 


but  sin  29-2  sin  9 cos  9 = 


cos  29  - cos  ^ 9 - sin^  9 - Y 


a.xY  _ a.XY 

X. 


and 

I 


r - (x\r')Jb„+(x-Y^)c72.f-2XYC^+  zXY('xVY")d2^. 


Now  obtain  r • 


^ - 2ttX  +2C.Y  + 

B X 


2b.+ + /2XYc/^. 


I * 

Next  it  is  desired  to  have  CTi  — — , • 

Iy  = 2b.Y- 2<7a.Y  + 2QX +(ZXh<ixYh2: 
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and 


/ Vc" 

0Y’'  ~ 2 •“  Z^^ 

For  the  problem  of  plane  stress  In  an  isotropic  mate- 
rial  the  following  relations  between  the  components  of  stress 
and  strain  hold: 

0 - L./'o' _ 

3Y  (27) 


1- 


:lu  _ J—  Txy 

~C 


To  obtain  u'  (x,y) 


for  the  disk,  integrate 


* 


dx  “■  El  ^ 

f . 

lx  ^ b.-2<3^-^ /2xr^»-'j[zt.+  za^T- /zxYc/J) 

which  becomes  after  rearrangement  of  terms 

-|-[2  h. 0 - V)  - 2 aji+v]  + I Z X Y T/  - V)J  . 


Hence,  by  integration 


(y; 


Similarly 


/ ^ 

It~e|  2b.i-2c(,+  i2xyd^-V(2b.-2aif  /2XYdj)l 


or 


:-^(2M'-'')  + 2cfy/#v)  1-  i2XYd//-v; 


+ v.Vx) 


Integrating 

V=^^[zM‘-^)+2ai(''+v)]Y  + 6xyV»6-v) 

If  Dossible  it  would  be  desirable  to  determine 
and  Vj,(x)  before  oroceeding  further. 

Using  the  tiilrd  of  equations  (2?)  gives 

^XY-  ijy^±  r' 


or 


G' 


Using  expressions  of  u'ano  v'  above  leads  to 


Txy~^  ^ d\/'  c(i-^» 

' E:  r 


^ -f 
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/ f 6(/-v^)(xVYMcy,.  ^ c/v/ 

^ V e 77  ^ j . 


Eq’jate  this  expression  of  ^xy  to  the  one  derived  from  the 
strdss  function. 

It  is  known  that 

r'  - - 

^X3Y* 

Differentiating  above  with  respect  to  y yields 

dX 

/ 


Equating  the  two  expressions  for  Txy  leads  to 

r'  C ^0-v>)(x*-+Y‘idz  j,.'^  . . _ 

U \ £ ■*■  37  + 4y°)  )4- 


This  requires  that  if  d^^  is  not  zero  that 


or  that 

G ('-V)  _ , 

e “ • 

Substituting  for  G gives 

EO-n^)  . 1-w  . 

20 E + ^ ‘ • 

Therefore  d^  must  equal  zero.  Furthermore  this  requires  that 

4^'  -f  d ul zc^ 
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or  that 


77 


Integrating 

^0— 

u.'-  KsY  + K,  . 

Again  no  Interest  is  found  in  a uniform  displaceir.ent 
so  and  are  set  equal  to  zero. 

Also  since  v must  be  an  even  function  of  x due  to 
symmetry  t}:en  must  equal  zero. 

Like'-lse  u must  ue  an  even  function  of  y due  to  sym- 
netj'y  so  k^  must  equal  zero.  Consequently 


^ ^ -+-  d Ut,  _ Q — _ 2.  Cx_ 

d X dy  G 

and  since  ^ is  not  infinite  then  must  equal  zero. 

•»hat  ai-e  the  exnressions  for  the  d Isolacements  u(x,y) 
and  v(x,y)  in  the  plate?  It  is  necessary  to  know  these 
disnlacements  in  order  that  our  last  two  boundary  conditions 
may  be  applied. 

Recalling  the  relations  (3)  between  the  stress  and 
strain  components  for  the  two-dimensional  problem  of  an 
orthotrooic  material  and  the  expressions  for  the  stresses 
( l6 ) , (17),  and  (18),  the  disnlacements  may  be  obtained 
bv  integration 

V '-J 
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3X  E.  ^ 


|^  = R(A( 


iVt 

E*  E. 


) w, 


(z,+w.^ 


b(^ 


2.  I 


^ \J 

+ E 


-f  ^ — 

t,/  Wi 


Ci-z  + Wj)  ^ 


E * 


By  integration 


^.  ^ + U.(Y)  . 

W-  V 


(28) 


Similarly 


I?  ^ 


■aY 


-^l^l  E,  /U.C£.^W.^) 

. D / ; +/56*''JyA/  — '\  1 Su 


After  inte£-ratlon 


= R(ft  ( ' -i:,,.,  j *6 


i-  V.(x). 


(29) 


The  third  equation  of  equations  (3)  is 

XY  , 

Taklnt^  the  required  partial  derivatives  of  u and  v one  has 


lx  Iy  - 


I 'Jtx\{  — I I tggVvyV 

ET„  lUsdw^U^i-wJ 


E,  Aw.Cz.-w.W  ^ . E,  Aw.CH,two/J 


d 4(x)  . d Uq(y)  _ _i_  T, 


d X 


- C 


»fT. 


JT'T 


This  can  be  rewritten  as 


/yx  _ r^r  (j»T  ^ flo(€.(*^€  T I -ho^  £*Vyx^ 

’xY- K^w,Cz,tvV.)v  El^  j 


■4l 


seE^i 


In  order  to  coiapare  this  expression  for  the  shearing 
stress  with  (l8)  it  is  convenient  to  change  the  forn  of 
this  expression.  i^ewrlteTxy  as 

/Y  _ of  Loief^Gg^ f 

\ E,  a<ie)^Erl 


-V 


(.  /3&6Gxy  / 6%\  V^rY  _ I f V 

WxfZi+Wx) 


r ^ dUo\ 


(The  reader  nay  wish  to  refer  to  pages  8 and  9 to  refresh 
his  memory  on  some  of  the  relations  vsed  below.) 


One  has 


Y'  +'^*Y  . i+o<eVYx^  — -L  * €.^\Ax\  r 

UA"Tr^ 


. ' Also 


'XT 


since 

Oyx  _ 
f-f  ■ 

so  that- the  expression  becomes 


or  Vv  - \) 




xY 
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Recall  that  o(^i3  — I 


dL\"E.\  ^ /3  ^ / 


Since 


ly  Ve.e"  /_1_  _ _ Ve>E<  _ \J«Y 

'^=S-^(g,.v  eJ-  2C„  ’ 


then 


= 2 K + = 


G 


»r 


^*^+5  + ZV«r 


or 


_ VEkEy 

e*- 


G 


»r 


'■Che  exoresslon  then  becomes 


G.y  >1^.  _ I \/K  _ ^ _ , 


6'Ey  G„ 


Sii.-.ilarly 


^ _ / f^^eVrA 

\ E,  Cc^6)"E,  I 


-I 


Therefore  written  much  more  simply  as 


^3 


Comparing  (30)  with  (l6)  it  is  apparent  that  the 
obtained  using  the  stress  and  strain  component  relation- 
ships is  equal  to  the  % 0 'derived  from  the  stress  func- 
tion if  and  only  if 


Since  G^ty  equal  to  zero  then 

zero  and 


dv4 


-f 


d 1^0 


equals 


dVo 

W 

dUo 

HT 


Therefore 

Vi=  K,X. 

and 

u,  - -K,Y  +«}  . 

Since  there  is  no  interest  in  the  case  of  uniform 
dixplacements  in  the  x and  the  y directions  let  k = 0. 

Consider  v;  it  is  obvious  from  symmetry  that  v must 
be  an  even  function  of  x,  that  it  is  impossible  for  v to 
be  a linear  function  of  x.  Therefore  k,  must  equal  zero. 

Likewise  u can  not  be  an  odd  function  of  y and  hence, 
if  it  were  not  true  already,  it  would  be  true  now  tliat  k, 
must  equal  zero. 


'■A 


It  has  been  shown  that  u^(y)  and  v^(x)  and  that  u^(y) 
and  v^(x)  are  identically  zero.  Also  and  d^  have  been 
shown  to  equal  zero. 

Setting  equal  to  u(x,y)  on  the  boundary  of  the 

disk  gives 

r/  \ "1  p r - ( > h I 

^ ~ ^ ^QOfo(t)Ccos6i  + (st«e')Ex 


_ (b,  I ^ S Q CO 

q(i*3€)(cos9  iLne)Ex\  E-x 


s e 


Recalling  that 


c^  = o ^ 


o<6  - 


_ - B- 


and  that  A 2_  ~ it  follows  that 


/ 

0(6  + I 


fiz  - 


^ ~ I 
/S9 


fl  B2-0, 


This  leads  to 


o<e-i 

0<t  V( 


TT  O 

(36  t I i ^ 


, . 0(6-1 

and  since  : — 

0(6  +- 1 

unless  o(  6 — 


_ 5^  -I 
ae  f I 

Q , then 


is  not  equal  to  zero 


Ai 


— O I consequently 

= O . 


As  a result  of  the  fact  that  - O and 

of  actually  taking  the  real  pai’t  of  the  expression  for  u it 
follows  that 

^(tbJi-^)-zaJit^))cose=  £.ose) 

^ ' ''  I 


y>cY  /-B.cose]  + ^ Q 

Q(l^3i.)E}  / E, 


Comparing  the  coefficients  of  cos  0 on  both  sides  of 
the  equation  gives 


/3\\  \). 


t 


(31) 


The  final  boundary  condition  requires  that  for  x = 
a cos  0,  y '=  a sin  0 that  v*=  v. 


0 B,  ^'rxYc.os  o-i  5(n0^  _ 5 V^x  <3  5/n0  ' 

E"^  j Fv 


It  follows  from  this  that 


■iS.f  + 01(1+1)) I r ^ _ J_L?lJihx_B  — 

^ ^ ' •<ca('i+«<€)E^'  a€Q(i-t3()Ey  ‘ 


(32) 
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It  is  now  possible  to  deternine  uniquely  the  constants 
A,  and  and  in  so  doing -to  completely  solve  the  prob- 

lem. 

Substitution  of  the  values  of  and  O 2^  from  equa- 

tions (2p)  and  (26)  into  equations  (31)  and  (32)  give  two 
equations  in  the  two  unknowns  and  . Since A, 
and  B = B,  , the  subscripts  will  no  longer  be  carried. 

Thus  /\  and  ^ can  be  uniquely  determined. 


+ 


(^^-0  A 


6S€~i)B 

Ua\3€ti)  H- j ^ 


^ 3 v(fY  g ^ 

a(it3€)E^ 


Rewrite  as 

f l~\J  I (<tV)(o<fe-|) 

\ZaE  2a(o(e-n)E  . Ua£  za(aeti)^  o(ita^)Ej 


\ 2E 


+ 


Qrii:i)i  _ s 

2E  eJ 


o. 


This  simplifies  somewhat  to 


V6--J 
a(«€f  i)E 


+ 


wVv-i).,  '\/y  Ji^>)  ■). 
aCi+o<e)Ej 


(33) 
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Substitution  into  equation  (32)  leads  to 


■f 


(^€-\)b 

tf-d^(6€  +/) 


■+ 


■ 


_ /-y-o(£U,  /\  __  g _ 5jy,g  ^ 


Regrouping  and  simplifying  leads  to 

( i-\)  (iiV)(o^€-i)  n-atgVy>  \a  , /ii^  .(i*we-i)  ^ i±ej\>^Q 

^20E  Z(3('o<£  + i)e;  za(6€i-i]E  sea(i*se)£^' 


ZE  2E 


- o 


or 


(-c^feV> 

a(i  + o{OE 


-f 


\ + ot'i  Vyx 
«v€a(i+o<6)E, 


j-^  + f a(-^  .'ii  i 

a(lf£e)E  /3€a{H-B€)Eyj  \ E.J 


Solution  of  equations  (33)  and  (34)  simultaneously 
yields 


A . -t-  (flASv,T)El]  - _ (o(e-n)(e-E,')Q*5 

(fl€ti  ^ -u)  E,  t (o(*£Vu„)E]  “ E,  +f«‘e'-»V„)  E 


l|/3 


* 0<  D _ ci<gOf.ie)(v£Y~^YxQ^  ^ 


■(E-E,)a"S  1 

'o^(]  +0<6)[(5e  -^€  Y^]Ey  t ( 1 £j 

{«€-V)E^^-ro^*6%V,y)E 

3(1  f Be)  po<6 -a^^v/)  + ( 1 1 c(V\)yg) Ej 

^ [g6(/->g6\^)EY ->-fitaVVYx)Ej 


^.,€Ot-^yv'e,-^)„E')  o'sl 

■(3<-v)E^+r-5’€'t^xY)E 

^f-o(*fV)E^  -t  (1  + «X*6’-\^yx)E 

(Be  m)  j?w  € -\>)  Ex  +C«<*^  V ^»y)  E] 

OL  [6€(i-S€\>)£,  (/3£'V)E^t(SW\),:;\£ 

a Cf  I '*<€  V)}  Ey+  (i + E]  (>ee 


q^£(\>Ey-\JYxE)  0^5  _ (E-~Ex)^^5 

_ __  o(e(>-«£v^)EY'*' (* (^^-'^)E)^’^ Co^^<£VMcv)E 
° ” ^ ZB€0-3€^)£^-^({+3*e\>rx)Z]_  C&e-^>)E^f(S"(£\\>^^E. 

6 r/+6e)[(«<^/ -c(ei;;EY-Ki+.(Vv|, Je]  I } [(«6 Ex-^f<><*e\-NJxY')E:] 


How  that  the  constants  liave  been  uniquely  determined 
a listing  will  oe  made  of  the  expressions  for  the  stresses 
and  the  displac ements  in  the  plate  and  in  the  disk. 

For  the  nlate: 


c<^£"A 

-A 

■+ 

W.CH.tW,') 

Wi(EitW 

(17) 


(16) 


h-9 


/SeB 


(18) 


u=r{a(-^^)(-^).b(- 


V^xy\ 

1 e:,  / 

\Zi+wJ 

(28) 


V = R ( , 8 (-^)(«fcsr) 


(29) 


_ SyU 


YVyx 


For  the  disk:. 

a;'=  zbo-^^i 

CTy  r zb.  + za^ 

XcT  — ® 

u'  --^(  X 

y*  ; + za^n+\;j}  Y 
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It  is  desirable  to  show  that  for  a rigid  disk,  that 
Is,  for  ET  =:  OO  , the  constants  A and  B reduce  to  the  con- 
stants obtained  by  Owens^^  for  a rigid  circular  disk  in  an 
infinite  orthotronic  elate  and  thus  that  tiie  solution  for 
this  condition  reduces  to  his  solution. 

Prom  (35)  one  obtains,  using  the  relation  cKS  r=  / , 

as  a limit  for  £ — >•  oo  . 

+ \ I Cl+«<‘eVyxXi5£t/X'^’^*^‘-  Xirt) 

5 ^ 1 1 ^e)(  I + o<'£  (I  tdf  ‘ V \^xt) 

Combining  the  denominator  into  a single  fraction  this 
becomes 

0-^5  a (i  +c(eX-^e\)yfi!'5)(Be\ v^;cy) 

Qc0t/g*£^\)yyyo<*6%  Vxy)  — 

/^Oi-/3eX I O^y) 


This  reduces  to 


€ Vyx')  + € V^YX  ( ‘ VxY^} 

<^0tx9"€Vy,)(o(*6%Mcy)--<5r5*e"+Vx>ryi-Ho<*'‘VYx') 


10 


Owens,  o£.  cit. , p.  32. 


which  is  the  value  for  A obtained  by  Owens  for  a circular 
disk  of  radius  a. 

Equation  (36)  gives  for  ^ ^ oo 


_ o<e^\  a^5 

/ V'0(*-£*'vlyv  o<*£*-t-V^; 


XT 


«<6+o<€)r/  f o-ho(€.)(/3^^^  v;,y) 


or 


B = 


^ /e<  € V>Yx  ( * ■»•«»< ’^xt)  -h  (|  » o<t)('l  VyJ 



<^(i**ieX  I \)^y)  -/9(^/fo^6X<gV%\jryYM-o<V\iYx^ 
;+/5€)('/fo<'-£‘VYxXo<*'e*'  + 


The  expression  for  B simplifies  to 

6 = -(K^£)o"sf— ■}. 

L ( 1 1 J(o<^£  *■+  vly)  Vw',y)0  ^C<  X,)  J 


This  is  the  same  value  for  B as  obtained  by  Owens  and 
thus  t'le  3ol\ition  found  above  oecomes  the  solution  obtained 
by  Owens  v/hen  E -DO  , It  has  thus  been  shown  that  Owens' 
problem  is  a special  case  of  the  above  problem,  namely  when- 
ever E - OO  . 

Furthermore  It  is  desired  to  show  that  the  problem  of 
this  dissertation  becomes  the  problem  solved  by  Professor 
Smith"^^.  It  is  recalle  i that  Professor  Smith  solved  the 


11 


Smith,  on.  c it . 


probleLi  or  an  infinite  ortliotropic  plate  ’.vith  an  elliptic 
or  circular  hole  at  the  center.  The  equality  of  the  tv/o 
solutions  will  be  shown  by  demonstrating  that  the  A and  B 
of  the  above  problem  become  the  A and  B of  Professor  Smith's 
problem  (a  = b)  whenever  E - 0. 

Equation  (35)  gives  for  E = 0 


{p<e-'0)  /$(i  f e 


(lt/96)Co(6-V) 


This  can  be  written  as 


!\  - -(v+e<€)0'5 


ocaeU-Be^))  ^ c^fcV(/36-V) 

(l->g£V)(o^£-V>)  — (^€-V)(|-o<£V) 

C/-o<ev')(«^€-V) 


Regrouping,  one  obtains 

A-  o<e-,8€-o(t\;*+s£v-) 

v/hich  reduces  to 


(\-\-aiQ.)g^5 
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This  is  precisely  the  A obtained  by  Professor  Smith  for  the 
free  circular  hole. 

Prom  equation  (36)  for  E = 0 one  has 


(H>crf€)  ^)Q^5  ^ (o(€  1 1)  Q^5 

I-  o^€y)  0(€.-V 

Oto(€)(i-^€V)  __  ^ I to(€)(a6-\>) 

atS6Xi-^ev)  0 + 


or 


gr  QfSe)a"'5 


\Jide  - VJ)  -t  1 - o(£v/ 

(I  - 66  V )(«< e 6 - VK  / -o<6  V) 


This  can  be  written  as 


B = (iiSe)a"5 


/-  

o(e  - -ci€V^-tS€^^ 


Therefore 

B = (i 

€ (o^--5) 

which  is  the  value  obtained  by  Professor  Smith  for  B.  Con- 
sequently, the  problem  of  this  dissertation  becomes  the  prob- 
lem treated  by  Professor  Smith  when  a = b and  when  E - 0. 


CHAPTER  III 


THE  ISOTROPIC  STRESS  FUNCTION  OBTAINED  AS  THE 
LIMITING  CASE  OP  THE  ORTHOTROPIC  STRESS  FUNCTION 


The  stress  function  for  an  infinite  isotropic  plate 
containing  a circular  Isotropic  disk  at  its  center  will  be 
obtained  as  the  limiting  case  of  the  orthotropic  stress 
function  (li].),  where  A and  B of  equation  (li;)  are  given  by- 
equations  (35)  and  (36)  respectively. 

It  is  a necessary,  but  not  a sufficient,  condition  for 
a material  to  be  Isotropic  that  E^  - Ey.  For  some  ortho- 
tropic materials  it  is  true  that  E^  3 Ey,  For  convenience 
it  is  assumed  that  this  equality  holds  here  since  this 
assumption  causes  no  loss  of  generality.  Thus 


= I 

a 


2. 


and  let 


C =; — i-uEit - £r_E 
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Q _ o(C  I ) 

(l  t Ey  + (» 

Q . e*^Cit-e”%v'E-VC') 

0 + C '^)CCe^-  e*»;')E  + ( 1 1 e^vi  E'] 

u _ (e-v)E,*(B^*\>„)e 

' (=t-  +(=<••+ Vx-,)E 

H - (g''^^-v')£  ^■  (e-‘^^.v)E' 

■ Ce‘^-v')E+  (c'»>+v)E' 

«(VEY-Vy,E)  _ e'^(v'E-VE') 

=^(i-«v)Ey  + (iv«*Vy,)E  e^(i-e'*iv')E  +(i  + e®V)E' 

1 E-E,  ^ C-E« 

' (e%.-  v')E  + (e^+ViE' 

|Y_  <xC('4lri-^V)Eyt(it3t«Y;)El  _ (^-v;E,  + (gVvj,y)£ 

/?( I t/S)[«(’l - X V)  Ey  + ( I + c**--a, El  (i  t V)  E,  + (xN  Vxy)E] 

I + e ®^)[e\i-e%')E -t (i+  e'^v)E’]  ^ H(e^V)Et(eW] 
f(2.)  :[;5:(Z,-W,)% 

<J(Zi,')  = E Ki  I n 


where 


2,  - X + I =<eY  r X + i.e%Y 

r?"  - K -V  i.  iSEY  rr  X H-  eY. 
^a.  “ 
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The  stress  function  P can  then  be  wriccen  as 


CD-GH 


f [h,c<^)] 


-h 


I- J 


cj[zS<^)\ 


■f 


Express  P as  follows 


5 pfO-e-»)(co-qH)  f[?^<»]  t(i-e* *)(i-j)q[?^(»)11[  ^ 5Y^ 
~ 2.  I (l-c*Xi-e"*')K  J ^ ' 


Por  an  isotropic  material  <P  becomes  equal  to  zero^^o 
As  (P~^0,  F becomes  indeterminate  and  assumes  the  form 


Let 

N = (l-e  VD-GH)ffc<'«lt(i-e®)fl-J)<3[?X‘»)] 

w=o-e*Xi-e'‘*“)K. 


Using  L'Hospital’s  rule  to  ascertain  the  value  of  P as 
^ O one  finds  that 

lim  N*  i lim  N' ' ; lim  W»  I lim  W* ' : 0 

<t>->o  0-»o  <?>-»o 

where  the  prime  designates  the  derivative  with  respect  to  <J). 

One  obtains  as  a limit  after  straightforward  but  some- 
what tedious  differentiation 

I/m  - _ 3 [(H-V')C  -t-(3-V)b1 

■ 0- v')E  + 0+v)E'  ■ 

^^Slncs  f : oo3h'’^5^  ( -i 

* ' ^XY  t-fc  J 


>"  (M 
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similarly  one  obtains  as  a limit 

lim  - _ :;^^r(2EVE)[0tV')E-Citv)E'] 

«’-*o  d 0^  \ +Ci+V)E']^ 

, CE-E)[(2t  v‘)E  -t-Cg-WE'T  -.■(VE-VE'^rv'Ei-f.-V^E'A. 
[(i-v)e+(hv)e']^ 

_ ^ [ E-  E -v'e  t VE']  / Ca'y  o*^  \ 

0-v»VE  + \ ^ 


Therefore  the  stress  function  for  the  isotropic  plate 
with  a circular  disk  of  a different  isotropic  material  at 
its  center  is 


‘feEVE)[^«»^*)E-0-»V)E*3-^CE-e)[(24-v')E»(s-V)E*] 
CCi-^')e:  +o^v)E'3[Cif\/')E  <-(3-V)E'3 


-h 


rv'E-VE')[vt  +(i-\/)E'] 

■ a ■ — 

[(■-V'JEf  Ci*V)E'JC(I+v'}E  <-(3-V)E'] 


2^itv)E-(i4V')E]  / u\  ^ ^ 

0+v;')E (3-v>;e'  V ^ J ^ 


-V- 
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APPENDIX 


For  a numerical  example  to  illustrate  the  theory  pre- 
viously developed,  an  infinite  Sitka  spruce  plate  with  a cir- 
cular copper  disk  at  its  center  subjected  to  a uniform  tension 
in  the  direction  of  the  grain  along  opposite  edges  of  the 
plate  will  be  considered.  The  X-axis  is  taken  in  the  direc- 
tion of  the  grain  and  the  Y-axis  is  taken  to  be  perpendicular 
to  the  grain.  The  normal  stress  components  in  the  plate 
along  the  two  axes  exterior  to  the  copper  disk,  and  the  shear- 
ing stress  components  along  the  boundary  of  the  disk  are  ob- 
tained. 

The  same  Sitka  spruce  plate  subjected  to  the  same  uni- 
form tension  as  above  but  with  a circular  rubber  disk  at  its 
center  is  considered.  The  same  stresses  as  above  are  obtained. 

For  an  example  of  an  infinite  Isotropic  plate  a steel 
plate  is  chosen.  It  also  is  subjected  to  a uniform  tension 
along  opposite  edges  as  previously  and  contains  a circular 
copper  disk  at  its  center.  The  same  plate  is  considered  with 
a circular  rubber  disk  at  its  center.  In  both  cases  the  same 
stresses  as  found  above  are  obtained. 

The  tables  and  the  graphs  are  self-explanatory. 
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TABLE  I 

ELASTIC  MODULI  OP  SITKA  SPRUCE"' 


E 


t 


G 


XY 


o,i|.6ij. 


1,679, 

000 

lbs. 

per 

sq. 

in 

76, 

000 

lbs. 

per 

sq. 

in 

112, 

000 

lbs. 

per 

sq. 

in 

The  elastic  moduli  of  Sitka  spruce  are  average  values 
obtained  from  a large  number  of  tests  made  at  the  U.  S. 
Forest  Products  Laboratory, 


TABLE  II 

CONSTANTS  EVALUATED  FOR  SITKA  SPRUCE 


^YX 

0.0210 

12.26 

€ 

2.168 

1.342 

4.700 

1.802 

K 

1.496 

/-f-dCG^Yx 

1.257 

06 

1.615 

4-  \)^Y 

12.72 

ex'*- 

2.608 

1 ^ a € v^Yx 

1.0378 

a 

0.6192 

3 ^ '^XY 

2.266 

0.3834 

1 + o<£ 

4.501 

oc€ 

3.501 

1 

2.342 

6o 


TABLE  III 

ELASTIC  CONSTANTS  FOR  COPPER,  HARD  RUBBER,  AND  STEEL 


E 

\) 

COPPER 

15,000,000  lbs.  per  sq.  in. 

Oo33 

HARD  RUBBER 

4,000  ' ' 1 t t t t 1 

0.25 

STEEL 

29,500, 000  ' ' 1 t 1 t ! 1 

0.29 

TABLE  IV 

CONSTANTS  A AND  B EVALUATED  FOR  AN  INFINITE 
PLAIN-SAWN  SITKA  SPRUCE  PLATE  CONTAINING  A 
CIRCULAR  ISOTROPIC  DISK  OF  RADIUS  a AT  ITS  CENTER 


A 

B 

COPPER  DISK 

0.3561  a^S 

- 0.1439  a^S 

HARD  RUBBER  DISK  - 2.0394 

1.0429  a^S 

TABLE  V 

VALUES  OF  THE  NORMAL 
X-AXIS  EXTERIOR  TO  A 
CENTER  AT  THE  ORIGIN 

STRESS  COMPONENTS  AT  POINTS  ALONG  WIE 
CIRCULAR  COPPER  DISK  OF  RADIUS  a WITH 
FOR  A PLAIN-SAWN  PLATE  OF  SITKA  SPRUCE 

X ; 

i 

a 1 

i 

1.1945 

0.0232 

2a  1 

1.1610 

0.0002 

3a  ' i 

1.1158 

- 0.0030 

4a 

1.0829 

- 0.0031 

5a 

1.0607 

- 0.0026 

6l 

TABLE  VI 

VALUES  OF  THE  NORMAL  STRESS  COMPONENTS  AT  POINTS  ALONG  THE 
Y-AXIS  EXTERIOR  TO  A CIRCULAR  COPPER  DISK  OP  RADIUS  a WITH 
CENTER  AT  THE  ORIGIN  FOR  A PLAIN-SA'WN  PLATE  OF  SITKA  SPRUCE 

Y 

O'vs 

a 

O.II4.O9 

0.0177 

1.05a 

0.531^2 

1.10a 

0.6812 

1.20a 

0.8091 

1.50a 

0.9214 

2a 

0.9656 

- 0.0065 

3a 

0.9869 

- 0.0029 

0.9930 

- 0.0016 

5a 

0.9956 

- 0.0010 

TABLE  VII 

VALUES  OF  THE 

SHEAR  STRESS  COMPONENTS 

AT  POINTS  ALONG  THE 

BOUNDARY  OF  A 

CIRCULAR  COPPFJi  DISK  OP 

RADIUS  a WITH  CENTER 

AT  THE  ORIGIN 

FOR  A PLAIN-SAWN  PLATE 

OP  SITKA  SPRUCE 

© in  de>a>ees 

^/S 

0 

0 

10 

- 

0.0007 

20 

000003 

50 

0.0050 

40 

0.0160 

50 

0.0372 

60 

0.0719 

70 

0.1162 

80 

0.1269 

85 

0.0831 

90 

0 
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TABLE  VIII 

VALUES  OF  THE  NORI.IAL  STRESS  COMPONENTS  AT  POINTS  ALONG  THE 
X-AXIS  EXTERIOR  TO  A CIRCULAR  HARD  RUBBER  DISK  OP  RADIUS  a 
WITH  CENTER  AT  THE  ORIGIN  FOR  A PLAIN-SAWN  PLATE  OP  SITKA  SPRUCE 


X 

a 

0.0115 

- 0.2024 

2a 

0.1210 

- 0.0252 

3a 

0.3574 

0.0061 

ij-a 

0.5373 

0.0110 

5a 

0.6600 

0.0104 

TABLE  IX 

VALUES  OP  THE 
AXIS  EXTERIOR 
CENTER  AT  THE 

NORMAL  STRESS  COMPONENTS  AT  POINTS  ALONG  THE  Y- 
TO  A CIRCULAR  HARD  RUBBER  DISK  OP  RADIUS  a WITH 
ORIGIN  FOR  A PLAIN-SAV/N  PLATE  OP  SITKA  SPRUCE 

Y 

_ . 05/s 

^/s 

a 
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Figure  Vsiriation  of  the  normal  stress 

component  at  points  along  the  X-axls  ex- 
terior to  a circular  copper  disk  of  radius  a 
with  center  at  the  origin  for  a plain-sawn  plate 
of  Sitka  spruce  and  for  an  isotropic  plate® 
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Figure  5. — Variation  of  the  nornial  stress 
component  at  points  along  the  X-axis  ex'^erior 
to  a circular  copper  disk  of  radius  a with  cen- 
ter at  the  origin  for  a plain-sawn  plate  of  Sit- 
ka spruce  and  for  an  isotropic  plate. 


66 


V 


Plgut^e  6o~  Variation  of  the  normal  stress 
component  at  points  along  the  Y-axis  ex- 

terior to  a circular  copper  disk  of  radius  a 
with  center  at  the  origin-  for  a plain-sawn  plate 
of  Sitka  spruce  and  for  an  isotropic  plate. 
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Flgiire  7. — Variation  of  the  normal  stresS' 
component  CT^  at  points  along  the  Y-axis  exterior 
to  a circular  copper  disk  of  radius  a with  cen- 
ter at  the  origin  for  a plain-sawn  plate  of  Sitka 
spruce  and  for  an  isotropic  plate. 
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Flgiore  3. — Variation  of  the  shear  stress  com- 
ponent points  along  the  boundary  of  a circu- 

lar copper^  disk  of  radius  a with  center  at  the  origin 
for  a plain-sawn  plate  of  Sitka  spruce  and  for  an 
isotropic  plate. 


Figure  9* — Variation  of  the  normal  stress 
component  at  points  along  the  X-axis  ex- 
terior to  a circular  rubber  disk  of  radius  a 
with  center  at  the  origin  for  a plain-sawn  plat 
of  Sitka  spruce  and  for  an  isotropic  plate. 
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Figur-e  10, — Variation  of  the  normal  stress 
component  0^  at  points  along  the  X-axis-  ex- 
terior to  a circulai’  rubber  disk  of  radius  a 
with  center  at  the  origin  for  a plain-sawn  plate 
of  Sitka  spruce  and  for  an  isotropic  plate. 
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Figure  11. — Variation  of  the  nornal  stress 
component  0^  at  points  along  the  Y-axis  ex- 
terior to  a circular  rubber  disk  of  radius  a 
with  center  at  the  origin  for  a plain-sawn  plate 
of  Sitka  spruce  and  for  an  isotropic  plate. 
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Figure  12. — Variation  of  the  normal  stress 
component  Gy  at  points  along  the  Y-axis  ex- 
terior to  a circular  rubber  disk  of  radius  a 
with  center  at  the  origin  for  a plain-sawn  plate 
of  Sitka  spruce  and  for  an  isotropic  plate. 
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Figure  13 . --Variation  of  the  shear  stress  com- 
ponent at  points  along  the  boundary  of  a circu- 
lar rubber  dis^  of  radius  a with  center  at  the  origin 
for  a plain-sawn  plate  of  Sitka  spruce  and  for  an 
isotropic  plate. 
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